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We study a one-dimensional Hamiltonian consisting of coupled SU(2) spin and orbital degrees of freedom. 
Using the density matrix renormalization group, we calculate the phase-diagram and the ground state 
correlation functions for this model. We find that, in addition to the ferromagnetic and power-law anti- 
ferromagnetic phases for spin and orbital degrees of freedom, this model has a gapless line extending from 
the ferromagnetic phase to the Bethe ansatz solvable SU(4) critical point, and a gapped phase with dou- 
bly degenerate ground states which form alternating spin and orbital singlets. The spin-gap and the order 
parameters are evaluated and the relevance to several recently discovered spin-gap materials is discussed. 
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In recent years a remarkable number of new materi- 
als have been synthesized, which despite having mag- 
netic ions have non-magnetic ground states, a finite cor- 
relation length and a gap in the spin-excitation spec- 
trum. Among them, the spin-ladder forming Stron- 
tium Cuprates H), the spin-Peierls systems Copper Ger- 
manates Q and the periodically depleted Calcium Vana- 
dates H have attracted the greatest attention. Here we 
discuss a novel spin-gap phase, which although math- 
ematically closely related to the spin-ladder and spin- 
Peierls systems, is physically quite different. It arises 
in systems with spin and orbital degeneracy and requires 
coherence over neighboring atoms of both degrees of free- 
dom. We also motivate such a spin-orbital model from 
a quantum-chemical analysis of recently discovered spin- 
gap materials Na 2 Ti 2 Sb 2 jj] and Na 2 V 2 5 §] and dis- 
cuss the relevance of this new phase to them. 

The two-band Hubbard model is well-known in con- 
text of magnetic insulators with Jahn- Teller ions |J . The 
Hamiltonian is: 

H = ^t°fcf a<7 c j [ )(T + ^2 UaprHourTHp*, (1) 
(a<7)/(/3<7') 

where i,j are site indices, a, (3 = 1,2 the orbitals and a, a 1 
the spin indices. Quarter filling of the bands amounts 
to one electron per atom. In the strong coupling limit 
this system is a Mott insulator, and the state of each 
ion can be characterized by a spin Si and the orbital 
state can be mapped into a pseudospin T — 1/2 so that 
orbital one corresponds to T z = 1/2 and orbital 2 to T z = 
— 1/2. Thus in the strong coupling limit, the effective 
spin-pscudospin Hamiltonian in one-dimension becomes 
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H — (JiSi • Si+i + J2T1 • T i+ i 

i 

+K{§ i -S i+1 ){f i -f i+1 )) (3) 

with Ji = J 2 = K/4. This special point has SU(4) sym- 
metry. However, there are many ways in which one can 
find deviations from these parameters. For example, it 
was shown by Arovas and Auerbach in the context 
of Cqo — TDAE, that more than one U a /3 leads to the 
Hamiltonian in Eq. 3 with the SU(2) x 5/7(2) symmetric 
parameter space reducing to J\ + J 2 = K/2, but with J\ 
not necessarily equal to J 2 . 

Here we study the ground states of this model with 
SU(2) x SU(2) symmetry with arbitrary J±, J 2 and 
positive K. The calculated ground-state phase dia- 
gram is shown in Fig.l, where the Arovas- Auerbach line, 
Ji + J 2 = K/2 is shown by the dotted line (which goes 
through point A). It is well known that the presence of 
orbital degrees of freedom can alter the nature of spin 
order, giving rise to both ferromagnetism and antiferro- 
magnetism |3). The ground states in phases I, II and 
III are known exactly. They are direct products of a 
spin and a pseudospin ground states. In phase I both 
spin and pseudospin degrees of freedom are fully polar- 
ized ferromagnets. In phase II the pseudospins are fer- 
romagnetic whereas the spins are antiferromagnetic and 
their ground state is the Bethe Ansatz ground state of 
the spin-half antiferromagnetic chain. In phase III, spins 
and pseudospins are interchanged with respect to phase 
II. All these 3 phases are conventional magnetic phases 
appropriate for ID, with gapless spin excitations. We 
will concentrate in the rest of the paper on the region 
IV, which we show has an excitation-gap, and the crit- 
ical line AB, where gapless excitations exist. There are 
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two other points in the phase diagram, where the ground 
state is known exactly. The point A has SU(4) symmetry 
and is also Bethe-ansatz integrable |8| . It has gapless ex- 
citations and power-law correlations. The exact ground 
state for point C was recently obtained by Kolezhuk and 
Mikeska Q . This ground state is doubly degenerate and 
there is a gap in the excitation spectrum. 

To determine the properties of the region IV and the 
line AB, we turn to numerical methods. We use the den- 
sity matrix renormalization group (DMRG) method 
to calculate the ground state, the excitation gap and the 
spin and orbital correlation functions. The overall fea- 
tures of the method remain essentially the same for this 
problem as for the Heisenberg spin problems Jll[. How- 
ever, to treat the biquadratic term properly, and to take 
advantage of the enhanced symmetry, modifications are 
needed. It is known that the product of operators within 
the same block in the block density matrix eigenvector 
basis will be error prone, because the resolution of iden- 
tity in the density-matrix basis is non- trivial ||. This 
creates a difficulty in rcnormalizing the biquadratic term. 
To circumvent this difficulty, we have carried through any 
operator products required to construct the Hamiltonian 
matrix at every iteration. Besides, as the system posseses 
a global SU(2) x SU(2) symmetry, we define a single site 
with 4-states from the spin and orbital degrees of free- 
dom and ensure that at every iteration both Sf ot and T t z ot 
are preserved as good quantum numbers. 

We have used periodic boundary conditions through- 
out this study and have verified our results extensively. 
The ground state energy for point C (see Fig.l) is cor- 
rect to numerical accuracy and the lowest excitation gap 
calculated to be 0.3756 in units of K compares with the 
variational estimate of 0.375 with the DMRG cut-off, 
m = 100. At the SU(4) point (point A in Fig.l), the 
ground state energy is accurate with that obtained by 
the Bethe ansatz [g| to fourth significant decimal place 
and the excitation gap vanishes to numerical accuracy. 
We also verify the ground state energies and quantum 
numbers of phases I, II and III and obtain their phase 
boundaries by DMRG. In all our calculations described 
below, we have kept the cut-off m to be 120 < m < 150. 
All the calculations are repeatedly checked by exact di- 
agonalization results for small system sizes and we report 
here the results for system sizes upto N > 30. Note that, 
each site contains a spin and an orbital variable. 

We have calculated the excited state energies in the 
subspace of (S 1 * 04 , Tj ot ) = (1,1) while the ground state 
remains in (S* ot , T 2 * ot ) = (0, 0) subspace. For the line BA 
as well as for the whole of the incommensurate line (line 
AD in Fig.l), we have calculated the excitation gaps from 
the N = An data. This is to avoid many-fold ground- 
state degeneracy for N = An + 2 Furthermore, to 
obtain these gaps in the thermodynamic limit, the calcu- 
lated finite-size gaps, for the N = An systems are fitted 
with a function of the form 



A(N) = A + A/N + B/N 2 + C/N 3 + .. . (4) 

We find that the SU(4) critical point A is the endpoint 
of a critical line AB, where the point B, is a special point 
where several ground states, including the fully ferromag- 
netic spin and pseudospin ground states become exactly 
degenerate. At the SU(4) point, we have verified that the 
spin structure factor peaks at q — n/2 and the decay of 
the power-law real-space correlations are consistent with 
the 3/2-power |14| . Furthermore, all along the open in- 
terval AB, the structure factor peak remains at q = tt/2. 
Also on this open interval, the ground state is a singlet 
and the spin and pseudospin gap remains zero. Numer- 
ically we find that the gap opens with an exponent of 
1.5 ± 0.25 along the line AD close to the point A. 

The point C is known to have degenerate ground states 
which can be written as a matrix-product consisting of 
alternating singlet bonds in both the spin and orbital 
variables |J . There is a finite correlation length and the 
ground state spin structure factor peaks at q — tt. Along 
the line AC above the SU(4) point A, the ground state 
remains doubly degenerate with a spin-gap (see Fig.l). 
This gap goes through a maximum at J = J\ = J2 = 
(0.5±0.02)i"T, while going from point C to point A where 
it vanishes to zero. Along this line, from point A to 
the point where the gap becomes maximum (point D in 
Fig.l), the spin as well as the pseudospin pair correlations 
are incommensurate. The peak in the structure factor 
moves from q = 7r at D to q = tt/2 at A. 

Consider now the line which runs from point C to in- 
finity with equal J\/K and J%/K values. All along this 
line, the spin-correlations stay peaked at (q — 7r) and 
there is a finite spin gap. In the small K/J limit, the 
results are similar to those found for the ordinary spin 
ladder, where a gap opens for any finite interchain cou- 
pling. Our results are consistent with those of Nersesyan 
and Tsvelik jlj| and Mostovoy in that the gap is lin- 
ear in K/J. In Fig. 2, we plot the lowest excitation gap 
in units of J — J± = J2, as a function of K/J. 

Over the entire region IV, the ground state with N = 
An is doubly degenerate. It is a spin and pseudospin 
singlet, with a finite excitation-gap. These results are in 
accordance with the Lieb-Schultz-Mattis theorem [p"7| . 
Furthermore, it has the same broken symmetry as the 
exactly solved point C. This is verified by calculating the 
square of the order parameter 

q 2 {i -j)= < (s?sf +1 - c)(s;s< +1 - c) 

+ (T?T? +1 -d)(T?T? +1 -d)> (5) 

for large i — j, where Sf and T? are the z-component 
of spin and pseudospin operators respectively, c and 
d are the average of two-particle correlations, i.e., c = 
UNEi < SfS? +1 > and d = 1/iVyV < TfTf +1 >, N 
is the number of spin-orbital pairs. We find that Q 2 re- 
mains finite in this phase. In the context of spin degrees 
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of freedom, it is appropriate to mention here that the 
critical line (BA) and gapful phases of this model dis- 
cussed above are similar to that of the Majumdar-Ghosh 
model for spin-1/2 chain with nearest and second neigh- 
bour exchanges Jl8[ . 

We now turn to some recently synthesized spin-gap 
materials for which this model is relevant. The material 
Na2Ti2Sb2 was recently found to have a finite temper- 
ature phase transition at T c w 1\QK Q|. Below this tem- 
perature the uniform susceptibility drops sharply with- 
out any formation of magnetic order, and shows acti- 
vated (spin-gap) behavior. In this material, the spin-half 
Ti 3+ (d 1 ) ions form a square lattice. Thus, this behavior 
is in marked contrast to the undoped Cuprates, which 
show antiferromagnetic order. 

This material has inverse "K2MF4" structure, consist- 
ing of layers of oxygen ions forming square-lattice, and 
Ti ions sitting in between them. The Sb ions are lo- 
cated above and below the centres of oxygen plaquettes 
as shown in Fig. 3a. Each Ti ion is surrounded by four 
Sb ions and two oxygens, forming approximate octahe- 
dron. In this local tetragonal coordination, the triply 
degenerate ti g levels are split into a doublet e g and an 
cti levels. As the covalency of the Ti— Sb bond gives the 
main contribution to the crystal-field splitting, the e g 
doublets should lie lower. Thus the ground state of Ti 
contains one d-electron in a doubly degenerate orbital, 
which choosing the O-Ti-0 axis as the z-axis are the 
d X z and the d yz orbitals. These orbitals overlap through 
strong 7T- hybridization with the corresponding p a orbitals 
of Sb thus giving rise to two independent one dimensional 
structures extending along two mutually perpendicular 
axes of the crystal (see fig. 3b). Thus to a first approx- 
imation this represents a quasi-lD quarter-filled 2-band 
system. 

Magnetic properties of the material Na2Ti2Sb20, are 
qualitatively consistent with the gapped phase discussed 
here, although, more experiments are clearly needed to 
investigate the relevant atomic orbitals [|l9| and the na- 
ture of the low temperature phase of this system. It 
is also interesting to note that the isostructural mate- 
rial Na2Ti2As20 (where As replaces Sb) in many re- 
spect shows similar behavior above a certain tempera- 
ture with a gradual reduction in magnetic susceptibility 
as the temperaure is lowered but also shows some evi- 
dence for ferromagnetism at very low temperature [p0| . 
This could be expected from changing parameters in the 
spin-pseudospin models discussed above. 

Another material that has recently gotten consider- 
able attention is Na2V2 0s. This material is also a two- 
band quarter-filled quasi-lD Hubbard system |21|. In 
this case, the two orbitals arise from two chains of a 
ladder-like structure. Hence pseudospin ordering corre- 
sponds to charge ordering on different atoms |22-p4]. It 



at the SU(4) point, the spin and pseudospin correlations 
are peaked at 7r/2 and this is different from the period 
1 and period 2 ( ferro and antiferromagnetic) charge or- 
dering scenarios, which have been discussed before. In 
the spin-gap phase discussed here, all atoms remain in an 
equivalent charge state. The question of whether such an 
alternating spin and orbital coherence between neighbors 
plays any role in this material deserves further attention. 

In conclusion, we have used the density matrix renor- 
malization group method to study the ground state, 
excitation-gap and correlation functions of an SU (2) x 
SU(2) spin-orbital model. We show that this model has 
a rich phase diagram. In addition to conventional ferro- 
magnetic and antiferromagnetic phases in the spin and 
pseudospin variables, there is a gapless critical line which 
runs from the fully polarized ferromagnetic phase to the 
Bethe-ansatz solvable SU(4) critical point. We also find 
that over a wide parameter range this model has a bro- 
ken symmetry gapped phase, where the system forms an 
alternating pattern of spin and orbital singlets. The rel- 
evance of this phase to various recently discovered spin- 
gap materials is discussed. 
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has been argued that in this material the spin-gap may 
be entirely due to charge ordering E4]. We note that 
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Figure Captions: 

Fig.l: Phase Diagram for the model in Hamiltonian (3) 
in the Ji/K, J2/K parameter space. Thick solid lines are 
lst-order phase boundaries and BA is a critical line. See 
text for details. 

Fig. 2: Lowest excitation gap in units of J as a function 
of K/ J for the line from point A to infinity through point 
D and C as in Fig.l. 

Fig. 3a: Structure of the material Na2Ti2Sb20. 

Fig. 3b: Cross-sectional view of the Titanium d xz and 
d yz orbitals, together with the p z -orbitals of Sb in the X- 
Y plane. As shown, hoppings t 11 (shaded-shaded)= t 22 
(unshaded-unshaded) = t and t 12 (shaded- unshaded) = 0. 
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